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Abstract 

We analyze the violations of linear fluctuation-dissipation theorem (FDT) in the coars- 
ening dynamics of the antiferromagnetic Ising model on percolation clusters in two di- 
mensions. The equilibrium magnetic response is shown to be non linear for magnetic 
fields of the order of the inverse square root of the number of sites. Two extreme regimes 
can be identified in the thermoremanent magnetization: (i) linear response and out-of- 
equilibrium relaxation for small waiting times (ii) non linear response and equilibrium 
relaxation for large waiting times. The function X(C) characterizing the deviations from 
linear FDT cross-overs from unity at short times to a finite positive value for longer times, 
with the same qualitative behavior whatever the waiting time. We show that the coars- 
ening dynamics on percolation clusters exhibits stronger long-term memory than usual 
euclidian coarsening. 
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1 Introduction 



Aging experiments in spin glasses [|], Q first carried out by Lundgren et al., have generated a large 
amount of experimental as well as theoretical work. Two types of experiments have been investigated: 
the zero-field-cooled experiment and the thermoremanent magnetization experiment, both leading to 
similar results. In the present article, we will restrict ourselves to the thermoremanent magnetization 
experiment, consisting in first quenching the system below its glass transition temperature at time 
t = 0, applying a small constant magnetic field up to the waiting time t w , switching off the magnetic 
field at time t w , and measuring the magnetization relaxation at time t w + r. It is an experimental 
observation that the magnetization relaxation depends on the "age" of the system, namely, on the 
waiting time. Different theoretical approaches have been developed so far, for instance: droplet 
picture [||, ||], mean field models || or phenomenological trap models ||. Several scenarios have 
been proposed, such as "true" versus "weak" ergodicity breaking Q, or "interrupted" aging ||. The 
first two scenario depend on whether ergodicity breaking occurs for finite or infinite waiting times. 
"Interrupted" aging means that, at a finite temperature, there is no more aging if the waiting time is 
larger than finite (but possibly large) time scale. In other words, the system equilibrates in a finite 
time. 

Aging can be characterized by the violation of the fluctuation-dissipation theorem (FDT). If the 
system has reached thermodynamic equilibrium before the magnetic field is switched off at time 
t w , the magnetization response is then independent on t w . This situation can be realized either at 
large temperatures, or in an "interrupted aging" situation (which will be the case in the present 
article), or, in the presence of non- interrupted aging, by formally taking the limit t w — > +00 before 
the thermodynamic limit N — > +00. Then, the equilibrium thermoremanent magnetization m(r) is 
related in a simple fashion to the autocorrelation of the spin configurations at times t w and t w + r 
via the FDT (see section || for more details). In the out-of-equilibrium dynamics, the FDT is no 
more valid, and there are analytical predictions in some mean-field solvable models of what is the 
FDT violation [7], ||. In particular, Cugliandolo and Kurchan ||] have proposed that the out- 
of-equilibrium linear response kernel R(t, t') relating the magnetization to the correlation depends 
on t and t' only through the autocorrelation C(t,t'). The FDT violation is then characterized by 
a function X(C) that depends only on the autocorrelation, and, as recalled in section ||, can be 
obtained from the thermoremanent magnetization simulations. 

The aim of the present article is to study the FDT violation in dilute Ising antiferromagnets at 
the percolation threshold, with a Hamiltonian 



where the summation is carried out over neighboring pairs of spins on a percolation cluster. In prac- 
tice, will will study here only percolation clusters generated on a square lattice in a two dimensional 




(1) 
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space. It is well-known that, for a finite cluster of N sites, the dynamics freezes as the temperature 
is decreased below the glass cross-over temperature || 

T - 2J2VP (2) 
±9 ~ IniV ' [Z) 

with d the fractal dimension and vp the percolation exponent. This glass cross-over originates from 
the conjugate effect of large-scale 'droplet' excitations (with zero temperature energy barriers that 



scale like Jin N [10, |llj), and the divergence of the correlation length at low temperatures JT^]. 

It is of interest to understand the FDT violation in these systems for two reasons. First, a quite 
different behavior from euclidian coarsening is expected, with more pronounced long-term memory 
effects due to the slow dynamics of 'droplet' excitations. We will indeed show that the function X{C) 
characterizing the fluctuation-dissipation ratio cross-overs from unity to a smaller value Xq in the 
aging regime. Whereas Xq is zero in euclidian coarsening, we find a non-zero value for coarsening on 
percolation clusters. This indicates that, even though these non-frustrated systems show interrupted 
aging, the FDT violation in these systems shares some common features with "true" spin-glasses. 
The second motivation for studying these systems is, as will be developed in section ||, that the low 
temperature magnetic response to an external magnetic field is linear only for magnetic fields smaller 
than a typical field h* scaling like T/yN (see section §). One is thus lead to study the FDT violation 
in the absence of linear response to an external magnetic field. 

This article is organized as follows: section [2] is devoted to analyzing the equilibrium magnetic 
response to an external magnetic field and to show that the low temperature equilibrium response 
is non linear. Section |3| recalls how the function X(C) characterizing the FDT violation can be 
obtained from the thermoremanent magnetization experiment. The results of our simulations are 
next presented and discussed in section ||. 

2 Absence of linear response at low temperature 

In this section, we analyze the low temperature equilibrium response to an external magnetic field. 
We consider a percolating cluster of N sites, and first analyze a toy model for the magnetization 
response to an external magnetic field. The equilibrium magnetization M(h) in an external magnetic 
field h can be expressed as 

d In (exp (/3/iM)) c 



M(h) = , ,3, 

where {X)q denotes the thermal average of the observable X with respect to the system without a 
magnetic field, and X denotes the disorder average. 



We are first going to formulate in section 2.1 a low-temperature toy-model which allows analytical 
calculations of (]3|). The predictions from this toy- model will be compared to simulations in section 
22, 
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2.1 Formulation of the low temperature toy-model 

Our toy-model relies on some assumptions about the geometry of the percolating clusters, and further 
assumptions about the low temperature magnetization distributions. The validity of these assump- 
tions relies on the fact that some of the predictions of our toy-model can be successfully compared 



to simulations (see section |2.3|) . 

In the dilute antiferromagnets model (|l]), the magnetization of the Neel state of the percolating 
cluster is equal, up to a sign, to the difference A = Na — Nb in the number of sites in the two 
sublattices A and B, the number of sites in the percolating cluster being N = Na + Nb- In order 
to allow for analytic treatment, we assume that both Na and Nb are independent variables and 
gaussian distributed according to 

with a width o scaling like a oc y/~N. Within these assumptions, the distribution of the Neel state 
magnetization A is also gaussian distributed: 

At zero temperature, the magnetization distribution of a given percolation cluster consists of two 
delta functions located at M = ±A. As shown in Fig. |], the effect of a small temperature is a 
broadening of the two peaks at ±A. The numerical calculations of P(M) shown in Fig. [l] were 
carried out using the Swendsen-Wang algorithm [13]. 

In our toy-model, we first make the approximation that all the geometry-dependence of the 
magnetization response is encoded in the single parameter A. This approximation becomes exact 
in the zero temperature limit. At a finite but sufficiently low temperature, we still assume this 
single-parameter description of geometric fluctuations. We further assume that the effect of a finite 
temperature is a gaussian broadening of the peaks at ±A in the magnetization distribution: 



The thermal broadening ap originates from low temperature excitations above the Neel state. At 
sufficiently low temperatures, only the lowest energy excitations contribute to a. If the system was 
not diluted, these excitations are properly described as a dilute gaz of clusters of spins with a wrong 
orientation with respect to the Neel state. The contribution to ap of these excitations is of the 
order of yNf(T), f being iV-independent and behaving like ln/(T) ~ —J/T. It is also well-known 
that long range low energy 'droplets' also exist in dilute percolating antiferromagnets, due to the 
self-similarity of the structure, and the fact that the order of ramification of the lattice is finite [14|: 
it costs a finite energy to isolate a cluster of arbitrary size from the rest of the structure. These 
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Figure 1: Magnetization distribution of a cluster of N = 836 sites, with A = 36. This cluster is shown 
in Fig. U (cluster A). 

'droplet' excitations can be clearly identified in the magnetization distribution of the ferromagnetic 



Ising model [15|. In fact, we can account for these droplet excitations in an effective distribution of 



the parameter ap over the Neel state magnetization A. We will come back on this point latter on in 



section 2.3 



We have chosen in (|6|) a gaussian contribution of thermal excitations. The resulting contribution 
to the magnetization of these thermal excitations is linear in the magnetic field. In fact, in order to 
describe how the magnetization saturates for magnetic fields scaling like TV , one should refine our 
toy-model to incorporate non gaussian tails in the magnetization distribution. Since, as mentionned 
previously, we are mainly interested in magnetic fields scaling like T/yN, these non gaussian tails 
do not play any significant role in this low magnetic field physics. 

2.2 Non linear effetcs 

Within this toy-model, it is straightforward to calculate the magnetic field dependence of the average 
magnetization for a fixed value of A. To do so, we notice that the magnetization distribution in our 
toy-model is nothing but the convolution of Pa,o and Po,ag ■ As a consequence, 

(exp(phM)) A ,* = (exp(/3/iM)) A ,o(exp(/5/iM)) , (T , 3 , 
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Figure 2: Variations of M{h)/a versus ha /T in the toy model calculation (see (§)), for various values 
of the parameter xp = Ja\/a. 



from what we deduce the average magnetization for a given value of A: 

x 2 h . , /hA 
M(/t) = ^- + Atanhf — 

For a fixed A, and for a magnetic field smaller than the cross-over magnetic field h* defined as 
the magnetization response is linear, whereas it is non linear for magnetic fields stronger than h*. 



Since A and ap scale like v iV, the cross-over field h* is small even for large systems. More precisely, 
we now average the magnetization over the geometry: 

/j r+oo / hA.\ h a f+°° /ha \ 

M(h) = a}- + J ^ P(A)Atanh [pjrjdb = aj- + ^= J ^ ntanh uj exp (-u 2 /4)du. (8) 

We should distinguish between the two regimes 

Weak fields: h < T/a M(h) ~ (2a 2 + aj)h/T 

Intermediate fields: T/a < h < Tcr/o| M(/i) ~ 2o-/^/vr + • 

As the magnetic field increased from zero, the response to the external field is first linear, and, for 
magnetic fields of the order of T/a, cross-overs to a non linear behavior. This behavior is shown in 
Fig. | for various values of the ratio xp = yo^/a. 
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Figure 3: Histogram of log ((h*) 2 A 2 /T 2 ), for different values of the Neel state magnetization A. The 
lines correspond to A = 1, 2, 3, 4 and the symbols to values of A between 5 and 19. The temperature 
is T = 0.4. 2000 clusters contained in the 20x20 square were generated for each value of A. 

2.3 Comparison with simulations 

We now compare the predictions of our toy-model for the equilibrium magnetic response of percolation 
clusters to numerical calculations. We have generated 2000 clusters for each value of the Neel state 
magnetization A. All these clusters are contained inside the 20x20 square. In order to compare with 
the toy- model results (^) , we have calculated for each of these clusters the cross-over field h* defined 
by the equality of the linear and cubic terms in the cumulant expansion (||): 



the magnetization distribution in a zero magnetic field being calculated with the Swendsen-Wang 



holds, at least for not too small values of A. This suggests that erg is not A-independent but is rather 
distributed, and scales on average like A. We attribute this behavior to the existence of low energy 
droplet excitations. Since these excitations correspond to large magnetic domains, the magnetization 
induced by reversing these domains should be coupled to the total Neel state magnetization A. 




0) 



algorithm |13[| . We have shown in Fig. || the histogram of log ((h*) 2 A. 2 /T 2 ). In the regime A S> 073, 
Eq. (^) becomes (h*/T) 2 = 3/A 2 : the histograms in Fig. |3| would be a 5 function located on the value 
log 3. Even if the histograms in Fig. ^ have a finite width, we see that the scaling (h* /T) 2 ~ 1/ A 2 still 
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3 Thermoremanent magnetization experiment 



In this section, we recall how the function X{C) characterizing the FDT violation can be obtained 
from the magnetization relaxation in the thermoremanent magnetization experiment. We first assume 
linear response. In the presence of a time-dependent magnetic field h(t), the linear contribution to 
the magnetization is 

m[h]{t) = I R(t,t')h(t')dt f , (10) 
with R the kernel response. We take here as a working hypothesis that 

R{t,t') =p X (C(t,t'))9(t-t')^fi, (11) 

a form of the response kernel suggested by the mean field model studies [||, 0, |8| . The auto-correlation 
C is 

C{t w + T,t w ) = (l^ C J i (^) C J i (t u , + T)). (12) 

In most of the spin-glass models (for instance: Sherrington-Kirkpatrick model, Edwards-Anderson 
model), the correlation length is vanishing, the reponse is purely local in space, and the terms i ^ j 
vanish in (|l^). In the presence of a finite correlation length, the thermoremanent magnetization 
is conjugate to the spatially non local autocorrelation (|i~2]). From this point of view, (10),(ll),(|l2]) 
can be safely taken as a extension of X(C) to our problem, in the sense that (i) lienar FDT reads 
X(C) = 1 (ii) we recover the usual definition of X(C) in the limit of a zero correlation length. 

Barrat used recently an interesting different method to handle spatially non local responses 
and non linearities: he measured the staggered magnetic response to a random field with a zero 
mean. In this way, the magnetic response to the random field is linear as a function of the width of 
the random field distribution, and conjugate to the local autocorrelations. We underline that Barrat 
does not consider the same conjugate quantities as ours, and the functions X(C) are thus different. 
In particular, we cannot handle symmetry breaking within our framework. However, in the case of 
the present problem of magnetic systems on percolation clusters, the non linearities are quite weak 
and, as we will see, we can characterize their effects on dynamics, which could not be possible in the 
framework of Barrat calculations. 

In the thermoremanent magnetization experiment, h(t) is a step function h(t) = h8(t w — t), so 
that the thermoremanent magnetization reads 

rC(t w +r,t w ) 
m{tw + t, t w ) = (3h I X(q)dq, 



where we have assumed C(t w + r, 0) =0 [we have indeed checked that this quantity was vanishing 
in our simulations]. The function X(C) is then obtained by differentiating the magnetic response 

T 

X{t w + r, t w ) = —m(t w + r, t w ) 
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with respect to the autocorrelation: X(C) = dx/dC. If the waiting time is large enough so that 
equilibrium has been reached, the magnetic response x(r) is i^-independent, X(C) is unity, and we 
recover the linear FDT: 

X (r) = ~m(r) = C eq {r). (13) 

Quite a lot of efforts have been devoted recently to characterize how the FDT is violated in an out- 
of-equilibrium situation. Analytical solutions were obtained in the framework of mean-field models 
H |7|, ||. The fluctuation-dissipation ratio was also obtained in numerical simulations in various 



models. For instance, in the case of spin glasses, Franz and Rieger [17] have analyzed the Edwards- 
Anderson model in three dimensions; more recently, Marinari et al. [18] have studied the FDT 
violation in three and four dimensional gaussian Ising spin glasses, and shown that the fluctuation- 
dissipation ration X{C) is, in these models, equal to the static Parisi function x(C). A model of 



fragile glass was also studied recently [19]. 



As explained in section [| the magnetic response of percolating dilute antiferromagnets is not 
linear for magnetic fields of the order of T j \fN . In the absence of linear response, the equilibrium 
magnetization can be expanded in powers of the magnetic field h, the coefficients of this expansion 
being the cumulants of the magnetization distribution (see Eq. (3)). Following the work of Gallavotti 



and Cohen pC|| , Kurchan [21] recently proposed an extension of the FDT to incoporate the effects of 
non linear response. However, we cannot use here this generalization in our Monte Carlo simulations 
since this would involve the calculation of the time-dependent magnetization distribution, included 
the tail where the magnetization is opposite to the magnetic field. For our purpose, we take here 
as a working phenomenological hypothesis that the thermoremanent magnetization is given by ([lC|), 
with the magnetic-field dependent response kernel 

,,dC(t,t>) 



R h (t,t') = [3X h (C(t,t'))0(t-t' 



Of 



In the presence of non linearities and out-of-equilibrium dynamics, Xh(C) contains contributions 
both from the non linearities and the aging dynamics. However, in the limit of large waiting times, 
the system has equilibrated and thus only the non linearities contribute to Xh(C). In the opposite 
limit of small waiting times, the non linearities do not contribute to X(C), and, in this limit, a 
contact can be made with FDT violations in other systems, especially euclidian coarsening dynamics 



4 Numerical results for X(C) 

We now present our numerical calculations of X(C). Our simulations were carried out on two clusters: 
a cluster A with N = 836 sites, A = 36, and a smaller cluster B with N = 294 sites and A = 16. 
These two clusters are shown in Fig. |j. 
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cluster A, 836 sites cluster B, 294 sites 



Figure 4: The two clusters that we studied. The cluster A contains N = 836 sites, and A = 36. The 
cluster B contains N = 294 sites, and A = 16. 



We first present in section |4,l| the equilibrium dynamics: the waiting time is long enough for 
the system to have equilibrated in the external magnetic field, and, on the basis of the arguments 
presented in section [2], we expect sensible non linear effects. In fact, the relaxation time is finite even 
in the thermodynamic limit (interrupted aging). As the size of the system increases, the relaxation 
time will first increase, due to zero energy barriers scaling like JlniV 0], and saturate when the 
linear size N l / d becomes larger than the correlation length given by [12| 



(2Jv P \ 
sT ~ exp [—jT- h 

with vp the percolation critical exponent. The limit of small waiting times is next presented in 
section [4.2[ In this situation, the magnetic response is linear. The combined effects of non linearities 
and out-of-equilibrium response arising for intermediate waiting times are next presented in section 
[4.3[ Finally, the dependence on r of the autocorrelation and the magnetic response are presented in 
section 



4.1 Large waiting times: equilibrium dynamics and non linear response 

We first examine the regime of a large waiting time t w , large enough for the magnetic response to be 
independent on t w . In practice, we systematically checked that the magnetic response was unchanged 
when the waiting time was increased by a factor of 10. The magnetic response Xh is plotted as a 
function of the autocorrelation C in Fig. || for the clusters A and B. We clearly observe on Fig. [5| 
important non linear effects since the magnetic response Xh depends explicitly on the magnetic field 
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Figure 5: Variation of the magnetic response Xh(C) versus the autocorrelation C in equilibrium 
relaxation. The temperature is T = 0.8, and the waiting time is t w = 10 5 . The magnetic fields are 
h = 0.1 (squares), h = 0.2 (circles) and h = 0.3 (crosses). The curves have been fitted along the 
procedure described in the text. 

h, even at the relatively high temperature T = 0.8. In the short time limit, we observe a behavior of 
the type Xh(C) = C — , whereas in the long time limit, Xh(C) = X^P'C. In order to interpolate 
between these two behaviors, we have fitted our numerical results to the form 

X h (C) = *|M = (1 _ xW) fc *, Xh (C) + X®, (14) 



with 

fc h x h {C) = (l + exp v ^ 

where \h controls the width of the cross-over between the short and the long time regimes, and 
C* = Cf } /(l - x£ 0) ). The fits obtained in this way are shown in Fig. [||, and, once the three 
parameters have been adjusted, a very good agreement with the simulation data is obtained. The 
variations of Xh(C) deduced from the fits are shown in Fig. || for the same simulations as in Fig. 
||. We observe in Fig. ^ that X^(C) cross-overs from unity at short times to a finite value X^ in 
the long time relaxation. If the response to the external magnetic field was linear, one would expect 
that X{C) = 1. Even though we could not address this question here, we expect a non linear FDT 
of the type pi] to hold in the long waiting time limit. 
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Figure 6: Variation of Xh(C) versus C in equilibrium relaxation. The temperature is T = 0.8. The 
variations of Xh(C) are deduced from the simulations presented in Fig. ||. 

4.2 Small waiting times: out-of-equilibrium dynamics and linear response 

In the short waiting time limit, the thermoremanent magnetization is linear as a function of the 
magnetic field. We have shown in Fig. |?j the variations of the magnetic response x versus the 
autocorrelation C for the two values of the mag netic field h = 0.2,0.3, and t w = 10 2 ,10 3 . Linear 
response is clearly observed. The variations of X(C) are shown in the insert of Fig. ^. Interestingly, 
the variations of X(C) in this situation where out-of-equilibrium effects are dominant are qualitatively 



the same as the ones in section 4.1 : X cross-overs from unity at short times to a finite value in the long 
time limit. We have no understanding of the reason why the variations of X(C) are qualitatively 
the same in the small and large waiting time limits, where deviations from linear FDT originate 
respectively from the out-of-equilibrium dynamics and non linear response. 

The fact that X(C) is finite in the aging regime is a quite noticeable difference with euclidian 



coarsening [16|, where X is vanishing in the aging regime (the linear response kernel R in (10) is 
zero in this regime). The long term memory of coarsening dynamics on percolating structures is thus 
stronger than for euclidian dynamics, which is maybe not surprising on the general grounds recalled 
in the introduction: a "droplet" of size N has a zero temperature energy barrier scaling like lniV 



, 10, 11], and a finite energy of the order of 2CJ, C being the order of ramification 
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0.0 0.1 0.2 



C 

Figure 7: Variations of the magnetic response xiP) versus the autocorrelation C in the short waiting 
time limit (t w = 10 2 , 10 3 ), and h = 0.2, 0.3. The corresponding variations of X(C) are shown in the 
insert. 



4.3 Intermediate waiting times: out-of-equilibrium dynamics and non linear re- 
sponse 

In order to examine the conjugate effects of nonlinearities in the magnetization response and out-of- 
equilibrium dynamics, we carried out the thermoremanent magnetization simulation with the cluster 
B at the temperature T = 0.55, and for a waiting time t w = 10 6 . The results are shown in Fig. || 
with r up to 10 7 . We have checked that equilibrium was not reached by carrying a simulation with 
a waiting time t w = 10 7 . On the other hand, the magnetic response Xh depends explicitly on the 
magnetic field h, as is visible in Fig. ||. We observe that X^C) can still be fitted by the form ([Ll|), 
even though we could not reach very small values of the correlation and magnetic responses, even 
for r = 10 7 . 

4.4 r-dependence of x(tw + r i tw) an d C(t w + r, t w ) 

In spin-glass models, the short time regime x(tw + T,t w ) = C(t w + r,t w ) — C^°'(t w ) is valid up to 
a time r* of the order of the waiting time t w fll~7|| . As shown in Fig. ^, we indeed observe such a 
dependence of r* in the out-of equilibrium situation: r* is of the order of 10 2 for t w = 10 2 , and of the 
order of 10 3 for t w = 10 3 . However, for larger waiting times, non linearities significantly reduce r* 



(r* ~ 10 3 for t w = 10 5 in Fig. ^ ). We observe in Fig. 10 in the case t w = 10 5 , T = 0.8 (equilibrium 
situation) that r* is also a function of the magnetic field h (the value of r* for h = 0.1 is one order 
of magnitude larger than for h = 0.2). This effect is also visible in the out-of-equilibrium simulation 
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Figure 8: Variations of the magnetic response Xh{C) versus the autocorrelation C for the cluster 
B, T = 0.55, t w = 10 6 , h = 0.1 (circles), h = 0.2 (squares) , 0.3 (crosses). The insert shows the 
corresponding variations of X(C). 
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Figure 9: Variations of x(^+r, t w ) +C^°\t w ) (diamonds: h = 0.2, crosses: h = 0.3) and C(t w +r, t w ) 
(solid lines) versus r. The waiting times are t w = 10 2 , 10 3 , 10 5 . The temperature is T = 0.8. 
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Figure 10: Variations of x(tw + T ^ w ) +C^(t w ) (squares: h = 0.1, diamonds: h = 0.2, crosses: 
h = 0.3) and C(t w + r, t w ) (solid lines) versus r. The waiting times are t w = 10 5 (T = 0.8), and 
t w = 10 6 (T = 0.55). 

shown in Fig. [10| (T = 0.55). However, from our simulations, we cannot make a precise statement 
on the variations of t* as a function of h for large waiting times. 



5 Conclusions 

We have thus carried out Monte Carlo simulations of the violation of the linear FDT in dilute 
percolating antiferromagnets. We have shown that these systems exhibit non linear response for 
magnetic fields of the order of T/y/N. In the small waiting time regime, the thermoremanent 
magnetization is linear in the magnetic field, but depends explicitly on the waiting time. On the 
other hand, for sufficiently large waiting times, the system has equilibrated (interrupted aging), and 
the magnetic response is non linear. Interestingly, in both situations, as well as in the intermediate 
situation where both out-of-equilibrium and non linear effects come into account, the function X[C) 
characterizing the deviations from linear FDT has qualitatively the same shape: it is unity at short 
times and cross-overs to a constant finite value for long times, the cross-over occurring at r*. In the 
small waiting time limit, r* is of the order of the waiting time t w . For larger waiting times, non 
linearities strongly reduce r* as the magnetic field is increased. By comparison with domain growth 
processes in non diluted lattices, the aging part of the dynamics shows stronger long-term memory, 
due to the existence of large scale low-energy 'droplet' excitations. 
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